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Abstract 

Random matrix ensembles with unitary invariance of measure (UE) are described 
in a unified way, using a combination of Tracy- Widom (TW) and Adler-Shiota-Van 
Moerbeke (ASvM) approaches to derivation of partial differential equations (PDE) for 
spectral gap probabilities. First, general 3-term recurrence relations for UE restricted to 
subsets of real line, or, in other words, for functions in the resolvent kernel, are obtained. 
Using them, simple universal relations between all TW dependent variables and one- 
dimensional Toda lattice r- functions are found. A universal system of PDE for UE is 
derived from previous relations, which leads also to a single independent PDEioi spectral 
gap probability of various UE. Thus, orthogonal function bases and Toda lattice are seen 
at the core of correspondence of different approaches. Moreover, Toda-AKNS system 
provides a common structure of PDE for unitary ensembles. Interestingly, this structure 
can be seen in two very different forms: one arises from orthogonal functions- Toda lattice 
considerations, while the other comes from Schlesinger equations for isomonodromic 
deformations and their relation with TW equations. The simple example of Gaussian 
matrices most neatly exposes this structure. 



*e-mail: igorrumanov@math.ucdavis.edu 
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I. INTRODUCTION 



Unitary ensembles of (single) random matrices (UE) are the most well known and thoroughly 
studied among random matrix ensembles (RME). Integrable properties possessed by their 
correlation functions and spectral gap probabilities, their ties with integrable hierarchies (see 
e.g. Ref. [10]) of partial differential equations (PDE) were extensively studied in the last two 
decades, see e.g. Ref. [23J. Still the approaches used to derive equations satisfied by gap 
probabilities, e.g. Refs. [29|[T|[8]. seem to be far from each other, moreover, the dissimilarities 
grow when considering more complicated RME like coupled or Pfaffian ensembles. Due to 
the aforementioned integrability, there are some PDE in every approach, which are universal, 
i.e. independent of the particular UE being studied. They may be equations satisfied by 
resolvent kernels of Fredholm integral operators and functions they are composed of. This 
was established as a general approach to study RME spectra in Ref. [29J (we call it TW 
further on). They also may be well-known nonlinear integrable PDE — the first members of 
KP or Toda lattice hierarchies, or their subhierarchies like KdV and others. This is the case 
in the approach started in Ref. [I] (called ASvM from now on), see also Refs. 0, El S ED] • But 
even this universal part of all methods, supplied by the integrable structure itself, and thus 
identical for different ensembles, is described differently in various approaches. As for the 
non-universal part, depending on the particular properties of a given probability measure, it 
practically leads to the necessity of almost case by case study of various ensembles. However, 
one might hope that it should be possible after all to describe the universal properties in a 
universal way, incorporating all existing approaches and emphasizing the common underlying 
structure inherent in all cases. 

First steps in this direction were made in Refs. [16] and [27], where the important TW 
and ASvM approaches were matched for several interesting cases: Airy, Bessel and Sine in 
Ref. [16J and Gaussian (or Hermite) in Ref. [27]. In the second paper some very simple 
relations between the dependent variables of TW and ASvM approaches for the Gaussian UE 
(GUE) were revealed. Here we extend these relations to all UE with spectrum on real line 
and expose some additional relations of the same kind, which seem to be overlooked in all 
previous considerations. (Everything we do here can be repeated for circular UE or even for 
orthogonal-polynomial ensembles with eigenvalues on a contour in complex plane, considered 
e.g. in Ref. [5J. Then, however, some specific formulas will differ from presented here, because 
of a different form of the basic three-term recurrence relations - our starting point, see below.) 
This allows us to derive some universal forms of PDE satisfied by gap probabilities of all UE 
irrespective of the specific spectral measure (or potential, in other words). Besides, we make 
contact with the isomonodromic approach [201 Ell EE EE El E] first applied to random matrices 
in Refs. [20~| Q31 [T3J H7J. In connection with TW work, it was most clearly exposed by Palmer 
in Ref. [23], where the generalized Schlesinger isomonodromic deformation equations for UE 
were written with a common form of non-universal, i.e. potential dependent, terms. We 
demonstrate that combinations of these last equations in general give the form of ASvM 
equations for r-functions obtained in Ref. [27] for the Gaussian case. These combinations 
expose an amazing similarity in structure with our universal PDE claimed above. At the core 
of all our considerations is still, as it was in Ref. [27] . the structure of orthogonal function 
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bases, three-term recurrence relations they satisfy, and one-dimensional Toda lattice hierarchy 
whose flows are in fact the continuous transformations among different such bases. We would 
like to stress here that while the last sentence is, strictly speaking, directly related only to 
finite size UE, the infinite ensembles can always be obtained as limits of finite ones. Moreover, 
the simple relations between TW and ASvM dependent variables, revealed here, allow one 
to define some quantities, originally meaningful only for finite matrices, like the r-ratios of 
matrix integrals of different matrix dimension, also for infinite ensembles. 

The plan of the paper is as follows. Section II gives a brief summary of the key results in the 
paper. In section III we derive the three-term recurrence relations for resolvent kernel related 
orthogonal functions. The functions are identified with orthogonal functions for ensembles 
with measure restricted to the complement of the original measure support [9] . Using these 
recurrences, in section IV we derive the universal relations between the TW and ASvM 
variables extending and generalizing results of Ref. [27J for the Gaussian case to all UE. 
In section V, from the results of previous sections, we derive new universal PDE satisfied 
by all UE. We obtain three equations reflecting the fundamental Toda-AKNS structure but 
in addition also another equation arising from the defining connection between logarithmic 
derivative of r-function and resolvent kernel. This allows us to obtain a single universal 
PDE for gap probabilities of UE modified by Toda times in ASvM approach, separately for 
each endpoint of the spectrum. Besides, we obtain a new interesting recursion structure for 
finite size matrix integrals. In the next section we introduce a different framework for UE - 
that of the isomonodromic deformations and Schlesinger equations, following the important 
work of Palmer [25J. We show then that our Toda-AKNS type system, derived by ASvM 
method in Ref. [27], can always be obtained also by summing up the Schlesinger equations. A 
partial summation of them gives also a Toda-AKNS-type system analogous to the one from 
previous section, this time, however, describing unmodified UE, i.e. ensembles with fixed 
couplings in the potential. The additional equation of previous section is now a consequence 
of summed Schlesinger system. Another consequence is two universal constraints for the non- 
universal terms in TW/Schlesinger equations. In section VII we consider simplest examples 
from the perspective of universal PDE (jl8|) of section V and explain the difficulty of its direct 
application to more complex cases. Section VIII is devoted to conclusions and open questions. 

II. BRIEF SUMMARY OF MAIN RESULTS 

For any UE on real line, given the corresponding orthogonal polynomial system and related 
Christoffel-Darboux kernel K n , functions Q n {x), P n {x) (n - matrix size), defining the resolvent 
kernel i? n of its restriction K£ to a subset J c of K also satisfy certain three-term recurrence 
relations, 



xQ 



n 



+ (a n + v n - v n+1 )Q n + 6 n _i(l - u n ){\ + w n )Q. 



n— 1 1 
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with coefficients explicitly determined by inner products u n , v n and w n ( 127j) over the subset, 
which were the additional auxiliary variables introduced by Tracy and Widom in Refs. [2E1 
to derive PDE for level spacing probabilities. Here u n = u n /b n _i, w n = w n /b n -i, and a n , b n 
are the coefficients in the original three-term relations (120]) . There are also simple recursions 
among the above quantities: 

Pn + l(x) = Qn-lfa) = (3) 

1 - U n l + W n 

W n +1 = U n /(1 ~ U n ). (4) 

Renormalized functions 

are orthonormal w.r.t. the original measure exp(— V(x))dx restricted to the complement J of 
the subset. Thus, they also satisfy a usual Christoffel-Darboux formula: 



R n (x,y) = f]r k (x)r k (y) = 6^(1 - fi^l + ^gi ^K-ify) r n ^(x)r n (y) 

and the three-term relations with symmetric (Jacobi) matrix of coefficients 

xr n {x) = b n (l - u n+1 ) 1,2 (l + w n+ i) 1,2 r n+1 (x) + (a n + v n - v n+ i)r n (x) + 

+ b n -!(l ~ U n ) 1/2 (l + Wn) l/2 r n -l(x). (7) 

The above formulas add in fact some more specific details to the earlier result of Borodin and 
Soshnikov [BJ on restricted RME. 

As a consequence of the above simple facts and general connection of orthogonal poly- 
nomial systems with Toda lattice, the simple universal relations hold with the UE matrix 
integrals, which are Toda r-functions: 

i rV) = 6.-1 ( i - r " +1 / r " +1 > ) = t r n t lTn : x ( i - r " +1 / r " +1 " 



U n = fa (I - K J n ) V) = K-l 



T nl T n J V ( T n) 2 V T nl T n 



Wn = (lP,(I-K n T )- 1 ^)=b n - 1 



T n-l/ T n-l 



T~n+l T n-l ( T n-l/ T n-l 



( T n) 2 V T n/ T n 



d t" 



v n = fa (I - J£)~V) ee fa (I - K J n )-fa) = -— In 



(9) 
(10) 



t=o 
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where t£ is the matrix integral over J n , r n - the same integral over R, and t± is the first Toda 
time. 

Due to the last relations, there is a universal system of PDE with respect to each single 
endpoint of the spectrum £ andt\ satisfied by T = lnr„ and r-ratios U n = and W n = 
if we consider the original measure modified by a linear potential, V(x) — > V(x) — t\x: 



— -f- = £— ? + 2-^C/„, 13 
d 2 W„ dW„ du 



*», ^ + 2 W W - <14) 
Here ([12]) is the well-known Toda equation, and the other three PDE satisfied by any such 
modified UE can be written equivalently as (t = t\, U — U n , W = W n below) 

d t {T^f = U ( Wx - WtU& (15) 
(T ?t ) 2 = -U ( Wt, (16) 

WU & - UW^t = £T ett . (17) 

It is possible to eliminate U and W from ffl5l) . ffl6|) and ffT2|) . and we obtain a universal PDE 
for the logarithm of gap probability P (T = lnr^ = InP + lnr n ): 

(T^ - r €tt T^ + 2(T Ci ) 3 ) 2 = (T c ) 2 (4T tt (T 54 ) 2 + (T^) 2 ) ■ (18) 

In addition, using also ( 117]) allows to derive a universal PDE for the logarithm of the ratio 
T + = U n = ln(r^ +1 /r I f) (where 'prime' stands for the derivative w.r.t. £): 

- - (n) a (r + )« - m ^" ]t + + (e - (r + )on(2(n)« - *) = °- ( i9 ) 

The universal Toda-AKNS-like structure of ( TT2]) - (|T4]) appears in a different guise from 
Schlesinger equations of isomonodromic approach, where Toda time t\ is not involved, in 
section VI. 
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III. RESOLVENT KERNEL\RESTRICTED ENSEMBLES' 3-TERM RELATIONS 

Consider a system of orthogonal functions (quasi-polynomials) {tp n } that satisfy a 3-term re- 
lation ( (p n (x) = TT n (x) exp(— V(x)/2), {ir n } are polynomials orthonormal w.r.t. a probability 
measure exp(— V(x))dx on real line R ): 

XLp n = b n cp n+1 + a n (f n + fe„_iv?„-i, (20) 

which can be written as 

xip n = (h<p) n , L = bA T + aA° + 6A, (21) 

i.e. L - infinite matrix with n-th row (L n ) 3 - = b n S n+ ij + a n 5 n j + & n -i^n-i,j (Ay = <5ij+i - 
the infinite shift matrix, a = diag(ao, di, 0,2, ■■•)> ^ = diag(bo, bi, 62, ■•■))■ The corresponding 
Christoffel-Darboux formula reads: 

^ / \ f{x)ij}{y) -ip(x)<p(y) ^ 

#n(z, y) = = > (p k (x)ip k (y), (22) 

y fc=0 

where <£>(a;) = ^/b n ^iip n (x), ip{x) = y/h^~Lip n -x(x) . Denote by K% the Fredholm operator 
with kernel fl22l) acting on £/ie complement J c of subset J of R. The resolvent kernel (we will 
sometimes use short-hand notation K for K^), R(x, y) = R^(x, y), the kernel of If (I — -ft') -1 , 
is defined so that the operator identity holds: 

(I + R)(I — K) = I. (23) 
Let us introduce the auxiliary functions playing a prominent role in TW approach [29J: 

Q(x; J) = (/ - K)-V(x), P{x- J) — (I — K)-V(x). (24) 

It is known since the seminal work [18] that if K is an integrable operator, i.e. its kernel 
can be expressed by the left-hand side of the Christoffel-Darboux formula (1221) . the resolvent 
operator is also integrable and its kernel can be written the same way in terms of the functions 
Q and P: 

p , x Q(x)P(y) - P(x)Q(y) Q n (x)P n (y) - P n (x)Q n (y) 

R n {x,y) = = 6 n _i , (25) 

x — y x — y 



where we let Q = y/b n -iQ n , P = y/h^P n , then 

<p n = (I- K)Q n , <^ n -i = (I - K)P n (26) 

Important for our purposes will be inner products u, v and w, introduced in Ref. [29] : 

u = u n = J JC Q(x; J)(p(x)dx, w = w n = J JC P(x; J)ip(x)dx, 
v = v n = f JC P(x; J)ip(x)dx = f JC Q(x; J)ip(x)dx 
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(27) 



the second equality in the definition of v being true by symmetry of the operator K£ and, 
consequently, also FL 3 n . Let u n = u n /b n _i, w n = w n /b n -i, then 

Lemma 1. 

n+l\ x ) — i =~ j tyn-l\ x 



1 - U n l + w n 

w n+ i = u n /{l-u n ). 
Proof. Functions P n+ i and Q n (as well as Q n -\ and P n ) are not independent: 

<p n = (I- K)Q n = (I- K n+1 )P n+1 = (I- K)P n+1 - (K n+1 - K)P n+1 , (21 
which can be rewritten as 



(J - (J - K)-\K n+l - K)) P n+1 = Q n . (29) 
+i 



On the other hand, the operator difference K 3 +1 — K 3 in (I2"9"j) is a simple projector (with = 



meaning the kernel of the operator on the left): 

K n+1 - K = (p n (x)(p n (y)xjc(y), (30) 

where Xj c (v) is the characteristic functon of J c , equal to 1 on J c and outside. After applying 
(I3"0"j) in (I2"9"j) and using the definition of quantities u n and w n , we find simple recursion formulas 
for P n+ i and Q n -i in terms of Q n and P n , i.e. ([3]). Shifting indices in fl3]), we find (jl]) as their 
consistency condition. 

□ 

Theorem 1. Functions Q n and P n in the resolvent kernel of a Fredholm operator associated 
with finite size UE, satisfy three term recurrence relations: 

+ {a n + v n - v n+1 )Q n + 6„_i(l - u n )(l + w n )Q n _ u (0) 

xP n = b n -i(l - u n )(l + w n )P n+1 + (a n _i - v n + v n -i)P n + 6 n _ 2 P n -i- © 
Proof. From ( 12 ip and (1261) we have 

xcp n = (L(J - K)Q) n = (I- K) (LQ) n - ([L, K]Q) n , (31) 
where the commutator gives 

([L,K}Q) n = (LKQ) n ~K(LQ) n = 
— {o n K n+ iQ n+ i + a n KQ n + b n _iK n _iQ n _i) — K(b n Q n+ i + a n Q n + 6 n _ 1 Q n „ 1 ) = 
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= b n {K n+1 - K)Q n+1 + b n _ x [K n ^ - K)Q n _ x . (32) 

Similarly, 



xif n ^ = (/ - KT)(LP) n - ([L, K]P) n (jMh) 

with 

([L, K\P) n = b n ^(K n+1 - K)P n+1 + b n . 2 {K n . x - K)P n - X . <m) 
Applying operator (/ — K)^ 1 to (l3Tj) . (l3Tk) we get 

(/ - K)-\ Vn = (LQ) n -(I- K)-\[L, K}Q) n , (33) 

(/ - K)- l x Vn ^ = (LP) n -(I- K)-\[L, K]P) n dMh) 

For the left-hand side of fl33|) we derive (imitating similar derivations of Refs. [HI [29] , see 
also Ref. [22]): 

(/ - KY x x^ n = -^=(1 - K)~ l x V = -^=(1 + R)x<p = 
1 ff s{x _ y) + Q(x)P(y)-Q(y)P(x) \ + {y _ dy = 



Jj\ x-y 
= (x{I + R)tp - J{Q{x)P{y) - Q{y)P{x))ip{y)dy\ = xQ n - vQ n + uP n , (34) 



\fK~- 

Completely analogously we get 

(I - K)~ 1 x(f n ^ 1 = xP n -wQ n + vP n . (35) 

On the right-hand side of fl33l) and fl33h). respectively, using fl32l) . fl32k). fl30|) and the defini- 
tions ff271) . we get 

(I - iC) -1 ([L, #]Q)„ = u n+1 Q n - b n ^u n ^P n , (36) 

(J - i^) _1 ([L, #]P) n = bn-tWn+tQn - V n ^P n . ([36b.) 

Finally, plugging ( |34|) . ( 136]) into (|33|) . and ( 135|) . (136k) into (133k). we come to the three-term 

relations ©, ©. □ 
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Consider now the difference of resolvent kernels for the operators of consecutive rank: 

R n+ i(x,y) - R n {x,y) = ^"j^^"^ = Q n {x)P n+ i(y) = P n+ i(x)Q n (y), 

l-u n 

therefore 

ry ( \ r Qn{x)P n {y) - P n (x)Q n (y) ^ 

Rn{x,y) = o n __i = } Q k {x)P k+1 {y). 

x — y 

y k=0 

If we introduce the following normalization for the orthogonal functions: 

/ \ Qk{x) ™. 

= (1-^)1/2 ' © 

then one can easily check that the new functions r k are orthonormal w.r.t. the measure 
exp(— V{x))dx restricted to J, i.e. the original measure on J and zero outside of it. Thus, in 
terms of these functions we can write out 

Theorem 2. The Christoffel-Darboux formula for the resolvent kernel holds: 

,1/2 r n(x)r n _ i 



n-1 

R n (x,y) = ^2r k (x)r k (y) = &„_i(l - u n ) 1/2 (l + w n Y 

k=0 y 



and the three-term relations with symmetric (Jacobi) matrix of coefficients are 

xr n (x) = b n (l - M n+ i) 1/2 (l + w n+1 ) 1/2 r n+1 (x) + (a n + v n - v n+1 )r n (x) + 

+& n -l(l - U n ) 1/2 (l + ^„) 1/2 r n _!(x). @ 

At this point one should recall the contents of the paper by Borodin and Soshnikov [9]. 
There the setup is completely identical to ours: two complementary subsets J and J c = M.\J 
of M are considered, and the corresponding operator defined by (3). Just the name "resolvent 
kernel" is not used, another name - "Janossy densities" - is applied. The authors, however, 
in fact give a simple proof of orthogonality of the above functions Q, P or r, so we just restate 
it briefly here in terms of Refs. [T8| 129] . By definition of functions Q n we can write 

n-1 

¥n = {I — K^)Q n = Qn~ ^k(^k, Qn)j<=, 

k=0 

or 

n-1 

Qn = ¥>n + 2_j Ank^k, (37) 
k=0 
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where A nk = ((fk,Q n )j c ~ the (matrix of) inner products, so that A nn = u n . Taking inner 
products of (1571) with ip k , k < n, over E, we see that 

(<fk,Qn)M = (<Pk,Qn)j<=, for k < 71, 

so {Q n } are orthogonal on J and 

(<Pn,Qn)m = (PnjPnjm = 1- 

It follows that 



(Qn, Qn) J — (<£>n, Qn) J — (y?n, Qn)]R ~ (<fn, Qn) J c — 1 ~ U n , 

which shows that functions {r n (x; J)} defined by (j5]) are indeed orthonormal on J C 1. 
Following Ref. [9], let us introduce the n x n matrices with elements (J, k = 0, 1, . . . , n — 1) 

(Gj) ifc = / ipj(x)ip k (x)dx, (Gjc) jk = 8 jk - (Gj) jk = / tpj(x)tp k (x)dx 
Taking inner products of (J57j) with over J c , one can easily see that 

8jk + A jk = (8j k - (Gjc) jk ) 1 = ((Gj)j k ) 1 , 

i.e. the relation for the restriction of to the n-dimensional subspace H = span{^}j = o,...,n-i 
equal to Gjc(I — G j<0~\ which was crucial in the argument of Ref. [9]. 
Remark. Also, see Proposition 2.5 of Ref. [7J, if we take the Gauss decomposition of Gj, 
Gj = S~ 1 (S T )~ 1 , where S is lower-triangular, then matrix S gives the coefficients of expansion 
of r k w.r.t. the (fj basis: r k (x; J) = X]j=o Skj{J)<-Pj{ x ), k — 0, . . . , n — 1. 

A more recent consideration of Janossy densities for general UE at the spectral edge, 
using Riemann-Hilbert problem techniques was undertaken in Ref. (26], also without making 
explicit connections with the framework of Ref. (29] . 



IV. UNIVERSAL RELATIONS BETWEEN TW VARIABLES AND TODA r- 
FUNCTIONS 

By definitions of the matrix integral over J n , t£, and the resolvent operator R 3 n = K^{I — 
K 3 ) -1 , we have the fundamental relation between them: 

RH(a j ,a j ) = (-iy- 1 ?^-. (38) 

This is just a consequence of the defining connection, the equality of two different expressions 
for the probability of all eigenvalues to lie in J, 

det(J - K J n ) = (39) 

Tn 

where r n is the corresponding matrix integral over whole M. n . Moreover, we have 
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Theorem 3. 



— I IT XfJ\-l \ L f -i T n+l/ T n+l\ jTn+lTn-1 ( \ T n+l/ T n+l 



Proof. Let us consider the values of our orthonormal functions r n (x) at the endpoints of J, 
cij. From the expressions preceding the definition (j5J) of functions r n and (1551) it follows that 

rn(«i) 2 = i^i(a ijai ) = (-ir 1 9111 ^ 1 ^ . (40) 

On the other hand, let us recall the universal TW equations [29] 

— = (-lYQ(a f ,jy, (41) 
^- = (-iyO(o J -;J)P(a J -;J). (43) 



From (141 D and (151) we have also 



r. W = (-ir' gj^jrJ . (44) 

Comparing the expressions (T4U]) and (|44l) . we get the simple universal relation (|S]) between 
w n and r-ratios. The other universal relation, for w n and corresponding r-ratios is obtained 
completely similarly, using (T42|) and Lemma 1. □ 

In fact, if we recall the recurrence relation between w n and u n -\ (J4]), we see that (J9j) is the 
same relation as flH]), just written in terms of w instead of u. The last two formulas were 
derived in Ref. [27] for the Gaussian matrices only, using the specific properties of this case. 
Now it is clear that they hold for all UE. 

For the Gaussian case we also obtained the following relation [27] : 

2v = 2(<p, (I - K J n )~^) = B-i lnr n J . (45) 

This relation, as we will see, is not true for general UE. However, it turns out that it has 
a universal analogue. To see it, the framework of Toda lattice for UE, established in ASvM 
approach [U [3], see also Ref. [30], sections 5-8, is crucial. Following these works, let us 
consider the UE with modified probability measure: p(x) = e~ v ( x > — > pt(x) = e~ v( - x * >+ ^ k = ltkX , 
t = (ti,t2,ts, ...), then the matrix integral 
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If n 
T n( t ) = -; A 2 n (x)T\p t (xi)dxi 
n\ Jjr, Y 



(46) 



is a r-function of integrable hierarchies - KP and 1- dimensional Toda lattice - as is the 
corresponding integral over whole M. n , r n {t). The orthogonal functions for t-deformed weight 
satisfy t-dependent 3-term recurrence relation 

xtp n (x; t) = b n (t)ip n+1 (x; t) + a n (t)ip n (x] t) + & n _i(i)</? n _i(x; t), 
and the coefficients are expressed in terms of the Toda r-functions: 



u\ d i r «+i( t ) 

an{t) = Trrln 7—, 

dti T n [t) 

T n +l(t)T n -i(t) 



(47) 



(48) 



Let us denote v + = v n+ \ — v n — a n , v _ = v n -i — v n + a n -i- Then 3-term relation (JTj) together 
with the same Toda lattice correspondence for restricted ( "resolvent" ) ensemble means: 



Tn+l^n—l 



->■ -M„)(l +W„) 



'ra+1 'n-1 



T. 



J^2 



g ln r "+l 



t=0 



-^ln r " +1 



-Am- 



T, 



n-1 



t=0 



(49) 
(50) 
(51) 



t=o 



dt 1 

From the last relations one can immediately conclude that 
Theorem 4. The inner product v n (t) can be expressed as 

and the universal relation for the original (unmodified) ensemble follows by taking the point 
t = 0: 



Otl T n 



(USD 



t=0 



Now, if we continue to consider the t-modified RME, the above universal relations together 
with 3-term recursion relations for the restricted measure (or resolvent kernel) ensembles will 
lead us to some universal PDE for matrix integral r-functions, containing derivatives w.r.t. 
the endpoints and the first Toda time derivatives. 
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V. UNIVERSAL PDE FROM 3-TERM RELATIONS FOR RESOLVENT KER- 
NEL 



Let us consider our 3-term relations from section 2 taken at a spectral endpoint £. Unlike the 
notations in Ref. [29] we denote q n = J), p n = P(£; J), i.e. the subscript refers to the 
size of the matrix rather than to the order number of the endpoint as in Ref. [2H]. Also let 
Qn = Qn(£,', J), Pn = Pn(£] J) ( see (ESJ), fl2"6"l) ). Then we have from 3-term relations ([T]), ([2]): 

b n q n+ i = (£ - a n + f„+l - W n )gn - &n-l(l - «n)Pn, (53) 



b n _ 2 Pn-l = (£ - Ctn-1 + - Un-l)Pn - &n-l(l + (54) 

Theorem 5. For all finite size UE, ratios of matrix integrals U n = t^ +1 /t^, W n = T^_ 1 /rj J l 
and function v n (t) = —dt\ ln(r^(t)/r n (t)) satisfy the following PDE w.r.t. any specific spectral 
endpoint £ G dJ and the first 1-Toda time t\: 

&Ur L _ dUn dVn ™ 

d*W n _ dW n dv 

Proof. Using formulas (jSJ) and (S3]), one gets for £ = a^: 

= (-l) fc 6 n gn + iPn + i = (-l) fc o n g n+1 -^-, (55) 

^2=1 = (-l) fc 6 n _ 2 g n _!p n _i = (-l) fc o n - 2 p»-i 1 f n _ , (56) 
oa k l + w n 

Substituting fl53]) and fl54"j) into fl55]) and fl56l) . respectively, we get 



(-l) fe 6 n _!9 n p„ + (f - a„ + v n+1 - t» n )(-l) fe — 



dVn+l _ 

a£ 1 — it 

<9w„ , dln(l-« n ) 

- (£ - a n + u n+1 - v n ) — , (57) 



<% n _i <% n <91n(l + w„) ,_ oN 

-gT- = — + - an-l + - ^n-l) ^ • (58) 

We introduce U n = t^ +1 /t^ and PV n = t^_i/t^. Then, with the help of relations (jHJ), © from 
previous section, fl57|) and f l58]) . respectively, can be written as (recall that v + = v n+ i — v n — a n , 
v- = v n -i -v n + a„_i) 
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= " 2 ^ iy - (60) 

Since, by f )52|) . u + = —9a lnf/ n , t>_ = — dn lnW n , the last two equations are equivalent to [T51 
and[H □ 

Thus, on the one hand, there always is a universal system of PDE valid for any UE and 
very similar in structure to the system derived in Ref. [27] for the Gaussian case using the 
first equations of 1-Toda lattice hierarchy [21 El [5U] - the coupled Toda-AKNS system (see, 
e.g., Ref. [21] )• It consists of the Toda equation itself: 



d 2 In ri dv n d 2 In t„ 



u n w n , (H2D 



dt 2 dti dt 2 

and the analog of AKNS derived above - the equations f[T5j) . (THI) . 

However, another important universal equation follows from matching with TW approach, 
if we recall formula 29 for the diagonal element of the resolvent kernel taken at an endpoint £, 

K{^Q=Pn-Qjr-<ln-Qr- (61) 

Theorem 6. With respect to any particular endpoint £ of the spectrum for a finite size UE, 
the following universal PDE holds for matrix integrals - r-functions: 



Proof. If we substitute (ISTl) into the square of (1551) . use TW formulas (14T|) . ( 142|) . and then 
formulas (J5J) and ()9j) from previous section, we arrive at equation (II ip . 

□ 

This will make it possible to derive something even more interesting - a universal PDE 
for only one dependent variable, the logarithm of the r-function r^. From the TW equations 
alone one can immediately get the simple bilinear relation: 

dv n \ 2 du n dw n 

(62) 



From now on we will denote T = lnr^ and write U, W instead of U n , W n and just t instead 
of t\\ also we will write just f^ for df/d^, ft for df /dt etc. Using the universal relations (J5J), 
(Q and (152|) with r-functions, (162|) can be rewritten as 



The last equation is also a consequence of the above (TT5j) . (Tl4j) . Thus, (TT5j) and (TT4j) can be 
alternatively rewritten as a different pair of equations - the equation (TIB"]) and their other 
independent combination: 
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The equation fill I) also can be transformed with the help of universal relations into another 
form: 

d^f = u ( Wx - WiU& (JED 

which will be used in what follows. It is convenient to introduce functions G and H, 

G = WU e - UWz, (63) 

H = WU t - UW t . (64) 

We now proceed to prove 

Theorem 7. The logarithm of gap probability P (T = \nr^ = lnP + lnr ny ) for random matrix 
unitary ensembles (UE) satisfies a universal PDE: 

(T^u - T&Pm + 2{T it ff = (T^ (4T„(T et ) 2 + (T^f) . m 

It would be interesting to find out whether this equation appeared before somewhere in the 
studies of nonlinear integrable PDE. So far the author is not aware of this. 

Proof. In terms of G ( I16p takes the form (using the Toda equation to eliminate U and W): 

G 2 = AT tt {T it ) 2 + (T itt ) 2 . (65) 

Now turn to the equation (TIB]) . One can express the combinations of derivatives of U and 
W in terms of T and G. Using the ^-derivative of equation ( fl2|) and the definition ( 163]) . we 

get 

2 = T itt + G, (66) 

ZUWz = T m - G. (67) 
Differentiating f )66|) and ( l6Tj) w.r.t. £ and using ([TE]) gives, respectively, 

2WU^ = 2{T^ t ) 2 + T m + G^ (68) 

2EW K = 2(T#) 2 + T Ktt -G { . (69) 
We now multiply (1671) by (168]) and subtract the product of (166]) and (169]) to get 

2UW(WtUx - E^W K ) = TsuGz - (2(T^) 2 + T m )G, 
or, using ([121) and ([TBI) . 
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T&Gz - (2(T ?4 ) 2 + T m )G = -2T tt d t (T ( ) 2 . 
Differentiating fl65|) w.r.t. £ one obtains 



(70) 



GG ? = T m (2(T & ) 2 + T Ktt ) + 4T«7V%. (71) 

After multiplying (!70!) by G and applying expressions f|65|) for G 2 and (|7T|) for GG^ on the 
left-hand side some of the terms and common factors cancel out and finally we arrive at 

T^G = TttTxu - TtuTw + 2(T^) 3 - (72) 

Thus, we have two independent expressions for the function G in terms of the derivatives 
of T - equations (!72|) and f|65|) . Substituting G from (!72|) into (I65I) we obtain ([18]) . □ 

Equation ( TT7|) was not used to obtain (fl8|) . It contains the ^-variable in the coefficient, 
while ( JTB"]) has the feature that its coefficients are constant. It turns out that the additional 
information from ([17]) can be used to obtain another universal PDE, this time for the logarithm 
of the ratio T + = ln(r^ +1 /r^) (now 'prime' stands for the derivative w.r.t. £): 

Theorem 8. The logarithm of the ratio of matrix integrals T + = h\(r^ +l /r^) for random 
matrix unitary ensembles (UE) satisfies a universal PDE: 



{T'l) tt - " (T[) 2 (T + ) tt - 1 + Z) +,t + + {i - (T + ) t )T' + (2(T' + ) t - 1) = 0. 



Proof. It is a tedious calculation. First one substitutes W expressed from Toda equation (Q2 
and W% from ffTB"]) into the ^-derivative of (fl~2]) . This gives 



T m = T tt T' + ?t 



n 

n 

Besides, we use (|T3|) . writing it as 



2T 5t = -{T' + ) t + (£ - (T + ) t )T|, 

and its t- derivative, 

2T 5tt = -{T' + ) tt + (£ - (T + ) t )(r;) t - T|(T + ) tt . 
We substitute the last two equations into the previous one, divide it by T+/4 and get: 

n 

We differentiate the last equation once again w.r.t. £ and compare the two expressions for 
T^ tt in terms of derivatives of T + . This gives ([TO]) . □ 



m t = -2^ - 2{T + ) tt + (r +)t y. 
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In addition to the main results - universal equations ( TT5|) - (IT7j) and especially (TT51) and 
(fH?|) . it seems worth to reformulate the preceding findings in this section. In terms of G and 
H and the main function T (iTTj) acquires a nice linear form: 

G t + Hz = 2^ tt . (73) 

Besides, if we use the three equations (IT2"l) - (fT4"j) . eliminating the mixed derivatives U^ t and 
W^ t by substituting their expressions from ([TBI and (JHj), respectively, into flT2"j) . differentiated 
twice, once w.r.t. £ and once w.r.t. t, we get (using also t>£ = — T^) 

T^t = £G + t^Wi + - 4T tt T#. 
This in turn can be rewritten as an expression of H in terms of G and derivatives of T: 

v otT 2T a T£ttt — TutT^tt — 8(T tt ) 2 T^ t 

H = 2£1 tt £ • (74) 

We thus have general universal expressions for logarithmic derivatives of Toda r-ratios T n ±i/r n 
in terms of derivatives of the main (n- level) r-function (due to explicit expressions (I65p . ( 172"|) 
for G and (JZU) for H): 

2^1nr/ = ^±^, 2^1nW = ^^, 

-t t* 1 tt 

r>o l rr T m -\-H T ttt — H 

2d t \nU = — , 2d t In W = — . 

J- tt J-u 

Therefore we can write out two pairs of general recursion relations for the derivatives of 
UE r-functions (restoring the subscript n for the matrix size) 



2d^T n+1 = 2d(T n + «9 5 In d 2 tt T n + 



2rp J 



2d ( T n = 2d^T n+l + «9 € In d 2 tt T n+l 



- 7 ' _ G n+ i 

dttTn+i 
and 



2d t T n+l = 2d t T n + d t In d 2 tt T n + 

Oft-' n 



2d t T n = 2d t T n+1 + d t lnd 2 t T n+1 



,2 r _ H n+ i 
r) 2 T - ' 



where G n = G, ff n = if at level (matrix size) n, and G n+ i, if n +i are given by the same 
formulas (172p . (174p . only in terms of lnr^ +1 rather than lnr^ . The first equations in the pairs 
follow from previous formulas for U, while the second ones - from the companion formulas 
for W with the shift n — 1 — >• n. 
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VI. ISOMONODROMIC DEFORMATIONS, SCHLESINGER EQUATIONS AND 
TODA-AKNS STRUCTURE 



The direct connection of TW equations for random matrix UE with general Schlesinger equa- 
tions - the compatibility conditions for systems of linear ordinary differential equations (ODE) 
related to their isomonodromic deformations - was first pointed out in the Ref. [T7]. There 
the matrices of the quadratic combinations of TW variables Q(ojt), P{a,k) (&fc is the k-th 
endpoint of the spectrum) were introduced, 

a _, uk-if Q(a k )P(a k ) -Q{a h f \ . , 

Ak ~ { ij V P (*kY -Q{a k )P{a k ) ) ' {7b) 

and shown to satisfy the Schlesinger equations. A little later an important work by Palmer [25] 
appeared, where the isomonodromic deformations for UE were studied in details. Palmer 
considered the corresponding "Cauchy-Riemann operator" d z problem similar to a Riemann- 
Hilbert problem (RHP) (for different or more general considerations of RHP in the context 
of matrix UE see e.g. an excellent book [11] as well as the original papers [TH [151 H2J El 113]). 
He studied the deformation of a fundamental matrix F(z), 

F(z) ' * * 



1p 1p 

of solutions to 2 x 2 linear system of ODE - called "differentiation formulas" in Ref. [2H] (for 
another, more recent, study of general matrix UE based on deformations of such differential 
systems see Ref. [5]): 

M(z) = 4- ( B [ Z \ ) , (76) 



m(z) V C(z) -A(z 

where m(z) is a polynomial in z, the functions A(z), B(z) and C(z) are determined [5J [29] . 
by the orthogonal functions <p, ip (see fl22|) ) and the corresponding potential V(x): 

f / \ / / \ V'(z)-V'(y) , Viz) 
J z-y 2 

2^, V'{z)-V'{y) 



B(z) = / <p\y) ■ ^dy, 

z-y 



V'{z)-V'{y) 



C(z) = / r(y) ■ — —dy. 

J z-y 

The matrix F(z), whose determinant must be a constant, is normalized by the condition 
det F{z) = 1. The deformed matrix Y(z; a) is uniquely determined by analyticity requirement 
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outside the real axis, the asymptotics near infinity, Y(z;a) ~ F(z),z — >■ oo, and matching 
conditions on a finite number of cuts [aj, ctj+i] on R: 

F + (^) = r_(2)-(/- X[aj , aj+l] + X[aj , aj+l] e), e=(J 2 ™ A ) , (77) 

where i + and y-(,z) are analytic in the upper and lower half-plane, respectively, and define 
Y(z) in their domains of analyticity. As shown in Ref. [25J, the solution to the above problem 
can be expressed by the formula 

y=d r x Si), P8) 



P ij-XKP 

where Q and P are the functions in evaluated for the rescaled operator \K J , i.e. Q = 
(I — Ai^ J ) _1 (/?, P = (I — XK J )~ 1 ip, and if is the operator with singular kernel (i.e. turning 
to 00 as x — > y) 

i> ( \ jOgMj/) - <p(x)il>(y) i x 

K[x,y) = XjW- 

x-y 

The normalization det F = 1 implies, by analyticity of det Y(z) and their asymptotic matching 
at infinity, the corresponding condition for Y, detY(z; a) = 1. This leads to a nontrivilal 
relation for the elements of Y: 

det Y = i)Q - <pP + A(P ■ KQ - Q ■ KP) = 1. (79) 
Then, following Palmer, introduce a matrix function A(z) such that, by definition, 

F{z)Y-\z) = I + A(z). (80) 

Differentiation of (180]) w.r.t. z combined again with analyticity considerations then yields the 
following ODE satisfied by Y: 

^f-Y- 1 = (/ + A)- X M(J + A) - (/ + A)- 1 ^, (81) 
dz dz 

where M = M(z) is given by (|7T>|) . The matrix A can be explicitly expressed in terms 

of related functions (p, ?/>, Q and P. It follows from f l75|) . ( 1T9"|) and the definitions flHUj) and 
(JTZD, that 



A(z;J) = ^^dyl ). (S2) 

J z y \ -^(y)P(y) if>(y)Q{y) 

Consideration of monodromies around the ends of the branch cuts aj for Y and the singular 
point at oo leads to the conclusion (23] that there exists a matrix function H(z; J) analytic 
everywhere in z-plane except for the point at oo and possible zeros of m(z) (see flT6]) . f )8T]) ) 
such that 
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with matrices Ak given by ( 175|) . In the simpler case, when m(z) = 1, the function Tl(z) is 
polynomial in z. One can see from general equations (1HTI) and (183|) that the principal part at 



p. p. (Y'Y- 1 ) = p. p. (II) = p. p. ((I + A) _1 M(J + A)). (84) 
Also in general m (z)U(z) is a matrix polynomial having no singularities besides oo so that 

p. p. (m(z)U(z)) = m(z)U(z). (85) 

This in principle determines matrix IT in general from M and A given by (ITS!) and (18"2"|) . 
The equations for the derivatives of Y w.r.t. the endpoints cij follow also from (IHUj) : 

—FY~ 1 d a YY~ 1 = d a A(z; J) 

or 

d^y- 1 = -(/ + A) -1 ^^; J) = 0(2 _1 ) as ^ ^ oo, 
due to the asymptotic expansion as z — > oo, 

A(z) = ^ + % + . . . (86) 

2 Z* 

The last equations determine that near the endpoints a,- 

p. p. (d.YY- 1 )- U ' 



2; — dj 



then by analyticity of d a YY 1 in z except for the poles at a,-, one has 



d a YY- 1 = (87) 



Z-dj 



Calculating the mixed derivatives d 2 za Y and d\ ka Y from equations and fl8T|) for different 
j and equating the singular parts when considering various limits z — >• a 5 -, gives Schlesinger 
equations modified [25] by II, see also Ref. [8]: 

^ = [n(a J ,J),^]-A^l^M, (88) 

*k = aI^I = (89) 
Another nonstandard Schlesinger equation for II arises from the same calculation [25] : 
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dU(z; J) _ [U(z; J) - IL( aj ; J), Aj] 
— A 



da -i 



z — a^, 



(90) 



3 * 3 

If we consider the limit z — > oo in the last equation as Palmer did and write M(z) 
M (z)(I + o(l)) near infinity, then, using f )84|) and the expansion f )86|) . we get 

[M ,^-A^]=0. 
Since matrix Mo is arbitrary, a completely universal equation follows: 

dAi 



with Ai determined by expansion of the formula 



Ai = A 



— v u 
—w V 



(91) 



(92) 



We see that parameter A drops out of (88) and we obtain exactly all three universal TW 
equations (141 p . f)42p and (1431) from Ref. [22] for the "auxiliary" variables u, v and w, which 
now are given the special importance and meaning in terms of r- functions, see (jHJ), ([9]) and 
([TO]) . There is also a general expression for logarithmic derivatives of w.r.t. the endpoints 
in terms of Aj and II in the isomonodromic approach (e.g. in Refs. [25| IS]): 



<91nT^ 
da* 



Tr 



a k 



(93) 



which is a consequence of existence of fundamental closed form Q, found first for the most 
general case of systems of linear ODE with rational coefficients in Ref. 



n = d a \nr J = d a \ndet(I - K J ) = J^Tr 



n(a i; J)Aj + 

k^j 



AkAj 
aj — at 



daj. 



(94) 



If we now sum up the equations (ISSj) and (l89|) over all endpoints aj and over all matrices 
Aj and use the equations ( 191~|) . we arrive at a system of three second order PDE for three 
functions u, v and w (the two equations for diagonal components are the same due to the 
tracelessness of the matrices involved): 



Bl 1 u = -^2\Il(a j ;J),A j \ i3 , 

j 



(95) 
(96) 



B\v = -Y J [n{a 3] J),A J } ll . 



(97) 
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Rewritten in terms of r-functions with the help of the universal relations (jSJ), ©, and (|52[) 
this becomes the universal analogue of the system of ASvM type obtained in Ref. [27] for the 
Gaussian case. Note, however, that the equation (1971) becomes in fact the third order PDE 
after using (152]) . For the Gaussian case, though, it turns out to be the total derivative of the 
first integral obtained in Ref. [29]. And the first integral itself is identified by (j45|) with our 
"boundary- Toda" equation [27J. 

Let us now sum up only equations (j88p and (1891) with derivatives w.r.t. a single particular 
endpoint a, only, using again (l9T|) and (1921) so that 



A ^ da k ' \ -w v 

k k x 

Again only terms involving II remain on the right-hand side. Expanding these commutators 
we write down 

Proposition 9. The three independent components of the Schlesinger equations associated 
with UE, summed over all for every matrix Aj, have the form of three-term PDE, 

A. B _ lU={no)i ^ + 2{n+)j ^., (99) 

d dw dv 

_ e _ ira ._ ( n ),_ + 2(n_),-, (loo) 

d „ ,„ s du ,„ . dw 

aSJ B - = (n -»^ + (n+) ^' (101) 

where we denoted IL = H(a,j] J) and matrix elements of U as 11+ = IIi2 ; II_ = II21, Ilo = 

n u - n 22 = 2n u . 

The last equality is again due to II being traceless as a consequence of tracelessness of the 
original matrix M(z) of differential relations. 

One can see a striking similarity in structure of the last equations and equations ({12]) - 
(TH|) from previous section. Only the current equation (llOip looks like a derivative of the 
Toda equation (I12p . which it is indeed in the special Gaussian case. The Gaussian case is 
now seen as a very special in the respect that there the general structural correspondence of 
the two systems holds literally term by term. This is not true for the other UE but still the 
Toda-AKNS-like structure and three-term character are present in both forms as if such exact 
correspondence existed at some deeper level. We believe this must be actually true although 
the real explanation is lacking yet. 

Let us now consider further the simplest (though important) case of only one endpoint 
£ = a,\. Then £>_i = and we rewrite the system (llOip . ( 199|) . (llOOp applying also the 
universal expressions of u and w in terms of r-ratios. Then, if we introduce 

n + = ^n + , n_ = ^n_, 
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our system reads: 



% = n ^ - 2n+^, (102) 

W a = -IW 5 + 2fl_^, (103) 

u K = -n_t/" { + n+w { . (104) 

We can add to this system equation ([531 . which for the single endpoint case can be written 

as 

r c = -n ^ + n_ M? - n + w 5 = -n ^ - n_t/^ - fi + w^. 

Here so far there are only ^-derivatives and no t- derivatives, however, the elements of IT, 
although known in principle from Palmer's considerations (see (1841) and (|85|) ). still do not 
have nice explicit expressions in general. Therefore, in general we obtain only some universal 
constraints for the non- universal H-quantities. The system (I102p — (I104p . if treated as a linear 
algebraic system w.r.t. n , n + and n_, is degenerate, and equation (TT6l) . (fg) 2 = —UgW^, is 
also a consequence of (I102p . (11031) like it was for (I13p . (1141) . The other independent combination 
of (11021) and (I103p together with (l§3"h) bring back our universal equation (JTJJ): 

WtU# - UzWk = 2v f (-n vf - n_c/"f - ti+w^) = 2^T e . 

Thus one can get only two rather than three universal constraints for Il-quantities, i.e. their 
expressions in terms of the derivatives of T = lnr;f. They are in fact equation ([93b ) and also 
(fTU4"P itself. 

VII. UNIVERSAL PDE ([18]) AND EXAMPLES 

Example 1: Gaussian emsemble. 

The Virasoro constraints [TJ, [3] give in this case (consider one interval situation - the distri- 
bution of largest eigenvalue): 

1 1 n 

dtT=--^T, %T = + 

We substitute these expressions for the derivatives of hir^ w.r.t. the first Toda time into our 
universal PDE (fT8|) and, denoting T' = d^T, get a fourth-order ODE: 

(T"T"" - (T'") 2 + 2(T") 3 ) 2 = 4(T') 2 ((T"') 2 + 4(T") 2 (T // + 2n)) . (105) 
This is actually a third order ODE for r = T'\ 

(r> r »> - {r»f + 2(r') 3 ) 2 = 4r 2 ({r"f + 4(r') 2 (r' + 2n)) . (106) 
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This third-order ODE is a nonlinear equation with coefficients not involving the independent 
variable £, unlike the third-order equation for Gaussian ensemble [2"9"] : 

r'" + 6(r') 2 + 8nr' = 4£(£r' - r), (107) 
which can be integrated once to give a form of Painleve IV equation [29] : 

(r") 2 + 4(r') 2 (r' + 2n) = 4(f r' - r) 2 . (108) 

One can verify that our equation (I106p is in fact a consequence of (jl08p . The Painleve equation 
(llOSp itself can also be obtained from (11061) . To this end define a function g{£) such that the 
right-hand side of (11061) equals 4r 2 g 2 and 

r'r'" - {r"f + 2(r') 3 = 2rg. (109) 
Differentiating the expression on the right-hand side of (11061) . defining g 2 (£) we have also 

(g 2 )' = 2gg' = 2r"r"' + 8r'r"(r' + 2n) + 4(r')V = 2r"{r"' + 6(r') 2 + 8nr'). (110) 
Eliminating the third derivative r'" from (11091) and (1 11 01) gives 

r 99 — r 9 + ^ rr 9 

or 

r' g ' - r "g = 2rr". (Ill) 

Substituting g = r'$ into (II 111) and solving the ODE for $, then using the boundary condition 
r — t- as £ — > oo, one gets 

g = 2(tr'-r), 

and, finally, using again the definition of g 2 from the right-hand side of (11061) . one arrives at 

(HDHD. 

Example 2: Laguerre ensemble. 

Here the potential is V(x) = x — alnx, and the relevant Virasoro constraints PQ, [3] give 
expressions of Toda-time derivatives in terms of the endpoint ones: 

8 t T = -£T' + n{n + a) , d tt T = £ 2 T" + n(n + a), 
so, for the function r = T', one finds 

r t = T (t = -tf - r, r tt = T m = ^r" + 2£r', 
r & = = -£r" - 2r', = T m = £V" + ^r" + 2r' . 
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For the Laguerre case it is convenient and standard to introduce the function a = £r instead 
of r, then one can see that 



/ /-a a /■ in i 

r t = -a , r tt = £cr , r (t = -a , r iu = f er + a 

and the universal PDE ffTBl reduces to the following ODE for o: 



e (-</(£</" + <0 + ^K) 2 - 2(^') 3 ) 2 = * 2 (W - o- + n(n + a))(a') 2 + £ 2 (a") 2 ) . (112) 
Again let us introduce function g(£) such that r.h.s. of (I112p = a 2 g 2 and 

Z (*W + °") - t^'f + 2(^) 3 ) = °9- (113) 
Differentiating the expression for g 2 on the r.h.s. of All 2j) gives 

gg > = " (£V + £a" + 6£(a') 2 - 4(7(7' + 4n(n + a)a') . (114) 
Again combining (I113p and fjl 14|) to eliminate the third derivative a'" leads to 

0*99' = °"g(g + v)> 

or 

a'g'-a"g = gg'\ (115) 
which is the same as fillip up to a factor of 2 on the right-hand side. So one gets 

g = io'-o. (116) 
Finally from the r.h.s. of (11121) and (11161) we arrive at 

eW') 2 + 4(£<t' - a + n(n + a))(o') 2 - (fr - a) 2 = 0, (117) 
which is the standard form of Painleve V equation for Laguerre ensemble [I] . 

The procedure just done for two classical ensembles can be applied for equation (ITS]) in 
general. If we put r = then the universal PDE (fT8"j) reads 

{nntt - runt + 2(n) 3 ) 2 = r 2 (AT tt (r t ) 2 + (r tt ) 2 ) . (118) 

We introduce function $ such that (actually r t § = G, with G defined by (1B"3"j) and satisfying 
(JSSD, (J72D) 

$ 2 = 2 

and 



$ = (f) + AT u (HQ) 
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r t r itt - ru^t + 2(r t ) 3 = rr t <5>. (120) 

Now differentiate flll9p with respect to £ and eliminate the most senior derivative r&t with 
the help of (I120p . Then one obtains 



*s = 77Tr < 121 > 



rr tt 

Two equations - (j!19p and (I12ip - together give an equivalent representation of the single 
equation (jll8p . If 012ip could be integrated in general to get an expression for $, then 
substituting it into (I119p would give a universal analog of Painleve equations for all unitary 
ensembles. Unfortunately, this is not the case. And, most important, there seems to be no not 
only general method to express the two i-derivatives - first and second - entering the universal 
PDE in terms of ^-derivatives and functions of £ only but even case by case application of 
Virasoro constraints beyond classical ensembles leads to infinitely many equations connecting 
derivatives w.r.t. all different times of Toda hierarchy, together with ^-derivatives. This 
seems to forbid simply eliminating the t-derivatives from the universal PDE for non-classical 
ensembles. To make this equation useful, a different point of view is needed. It might require 
some considerations of differential systems of the previous section, see also e.g. their treatment 
in Ref. [B], or a completely different perspective yet to be found. This open problem seems 
worth futher efforts since the approach presented here seems rather appealing. 



VIII. CONCLUSIONS 



The simple universal connections between eigenvalue spacing probabilities and their ratios 
represented as ratios of one-dimensional Toda r-functions on the one side, and auxiliary vari- 
ables naturally arising in the analytic approach to PDE for the probabilities using Fredholm 
determinants and resolvent operators on the other side, are derived. They hold for all Hermi- 
tian random matrix ensembles with unitary symmetry of the probability measure. They are 
obtained as a consequence of also derived here universal three-term recurrence relations for 
systems of functions orthogonal on subsets of real line. These orthogonal bases for restricted 
ensembles were studied by Borodin and Soshnikov in Ref. [9], so we make a missing connection 
of their work with Ref. [29J. All our considerations can be generalized to unitarily invariant 
random matrix ensembles with eigenvalues on a circle or other simple contours in complex 
plane. 

The above relations of functional-analytic approach with the structure of orthogonal func- 
tions and one-dimensional Toda (or Toda-AKNS) integrable hierarchy allow us to derive some 
universal PDE for 1-Toda r-functions-matrix integrals and their ratios. Even a single univer- 
sal PDE for the logarithm of gap probability (fl8|) of various unitary ensembles is obtained. 
Although the difficulty of applying the universal PDE f ll8p directly to obtain specific ODE 
for various non-classical ensembles is rather disappointing, the author believes that equation 
ffT8|) . and also f JT9|) for the ratio of consecutive size matrix integrals, may be important and 
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could be somehow used to describe both universal and model-specific properties of random 
matrix UE. How to use them remains an interesting question for future research. 
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